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Abstract 

The focusing of the vacuum modes of a quantized field by a parabolic mir- 
ror is investigated. We use a geometric optics approximation to calculate the 
energy density and mean squared field averages for scalar and electromagnetic 
fields near the focus. We find that these quantities grow as an inverse power 
of the distance to the focus. There is an attractive Casimir-Polder force on an 
atom which will draw it into the focus. Some estimates of the magnitude of 
the effects of this focusing indicate that it may be observable. 

PACS categories: 03.70. +k, 34.20.Cf, 12.20.Ds, 04.62. +v. 

1 Introduction 

The Casimir effect can be viewed as the reflection of vacuum fluctuations by mirrors. 
The presence of a reflecting boundary alters the modes of a quantized field, and results 
in shifts in the vacuum expectation values of quantities quadratic in the field, such as 
the energy density. Typically, Casimir effects for massless fields may be estimated by 
dimensional analysis. If r is the distance to the nearest boundary, then the Casimir 
energy density is typically of order r -4 times a dimensionless constant. This constant 
is usually of order 10~ 3 in four-dimensional spacetime. It is of course possible to find 
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a much smaller result due to special cancellations. For example, the Casimir energy 
density for a single, perfectly conducting plate is zero, even though the mean squared 
electric and mean squared magnetic fields are separately nonzero. 

These typical results arise from calculations of specific geometries, not from any 
general theorem. This leaves the possibility of exceptions, where the energy density 
is much larger than would be expected on dimensional grounds. Indeed, one possible 
mechanism for amplification of vacuum fluctuations has already been proposed |TJ : . 
This mechanism is based on the fact that the contribution of various parts of the 
frequency spectrum to the Casimir effect is a highly oscillatory function || Q]. The 
contributions of different ranges of frequency almost, but not quite, completely cancel 
one another. The possibility that one can enhance the magnitude of the effect by 
altering the reflectivity of the boundary in selected frequency ranges was discussed 
in Refs. |, |. 

However, in this paper, we wish to propose a different mechanism for amplification 
of vacuum fluctuations. This is the use of parabolic mirrors to create anomalously 
large effects near the mirror's focus. It is well known in classical physics that a 
parabolic mirror can focus incident rays which are parallel to the mirror's axis. This 
means that a particular plane wave mode becomes singular at the focus. One might 
wonder whether this classical focusing effect of modes can produce large vacuum 
fluctuations near the focus. We will argue that the answer to this question is yes. 

The outline of this paper is as follows: In Sect. [| the essential formalism needed 
to compute mean squared field averages in the geometric optics approximation will 
be developed. It will be argued that the dominant contributions will come from 
interference terms between different reflected rays. In particular, expressions will be 
given for the case of two reflected rays from a single incident ray in terms of the path 
length difference of the two reflected rays. In Sect. the specific case of parabolic 
mirrors will be studied, and the condition for the existence of multiply reflected rays 
given. It will be shown here that there is a minimum size required for a parabolic 
mirror to produce large vacuum fluctuation focusing. Section £| deals with a couple 
of technical issues, including the treatment of the apparently singular integrals which 
arise. In Sect. |5], we give explicit results in an approximation in which the mirror is 
only slightly larger than the minimum size needed to focus vacuum fluctuations. The 
possible experimental tests of these results are discussed in Sect. ||, and conclusions 
are given in Sect. ^. 

Units in which % — c = 1 will be used throughout this paper. Electromagnetic 
quantites will be in Lorentz-Heaviside units. 

2 Basic Formalism 

The approach which will be adopted in this paper is a geometric optics approxima- 
tion. This approximation assumes that the dominant contribution to the quantities 
which we calculate comes from modes whose wavelengths are short compared to the 
goemetric length scales of the system. The justification of the approximation will lie 
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in a self-consistent calculation leading to large contributions from short wavelength 
modes. At first sight, it might seem that this approximation would always fail, and 
that only modes whose wavelengths are of the order of the goemetric length scales 
will contribute significantly to quantities such Casimir energy densities. However, 
there is a circumstance in which this intuition can fail. This is when there are two or 
more reflected rays produced by the same incident beam. It then becomes possible to 
have an anomalously large interference term between these rays, as will be illustrated 
below. 

Let us first consider the case of a massless scalar field, (p. Let the field operator 
be expanded in term of normal modes as 

+ (i) 

k 

where a k and are creation and annihilation operators, and are the mode func- 
tions. The formal vacuum expectation value of ip 2 becomes 

<^ 2 >/ = EI^I 2 - (2) 

k 

In the absence of a boundary, the modes are simply plane waves. In the presence 
of the boundary, there are both incident and possibly one or more reflected waves for 
each wave vector k. Write the mode function as 

Fk = /k + E/k } , (3) 

i 

where /k is the incident wave and the are the reflected waves. (Note that here k 
denotes the incident wavevector.) We may take all of these waves to be plane waves 
with box normalization in a volume V, in which case 

^Tsr"^ (4) 

The take the same form, but with k replaced by the appropriate wavevector for 
the reflected wave. 

If we now insert Eq. (^) into Eq. (0), we obtain a sum involving both the absolute 
squares of the incident and the reflected waves, and the various possible cross terms 
between the different waves: 



(^)f = E 



IM 2 + E lif I 2 + + hf%) + E fPfK 



(5) 



This quantity is divergent and needs to be renormalized by subtraction of the corre- 
sponding quantity in the absence of boundaries. We will argue in Sect. 3~T that this 
is given by the above sum without the cross terms: 

(^ 2 )o = E(IM 2 + Eltf ) l 2 )- (6) 
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The renormalized expectation value is then given by the sum of cross terms 



(v 2 ) = (v 2 )f - (^ 2 )o = E 



E(/i/£° + /k/ w ;) + £/£W 



(7) 



The renormalization which we employ is the usual one of defining the renormalized 
quantity to be the difference between the formal expectation value with the mirror 
and that without it. Another way of expressing the same prescription is to say that 
quantities such as (ip 2 ) are only defined up to an additive constant, and we choose 
the constant so that (<£> 2 ) — > at infinite distances from the mirror. 
Let us examine a particular cross term: 
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E(tfV (2) k +/f/ (i) ; 



Here and are both of the form of Eq. ([|), except with k replaced by kx 
and k 2 , respectively. These might be two reflected waves, both corresponding to the 
same incident wavevector k, but different reflected wavevectors, ki and k 2 . Only the 
direction, not the frequency changes upon reflection, so 



|k x | = |k 2 | = |k| = u. (9) 



We can now write 



Tn = 2 Re £ * e <[-,- fa » _ _L / ,p k Mto. - fa) ■ xl (10) 

^ ZujV 87T J uj 

where the infinite volume limit has been taken. The argument of the cosine function 
is proportional to the difference in optical path lengths of the two rays, A£, so that 

(ki -k 2 ) -x = uA£, (11) 

and hence 

Tn = XlA^i. (12) 



87T 2 J UJ 

Note that the integral in the above expression will diverge as (A£)~ 2 in the limit 
that A£ — > 0. Thus within the geometric optics approximation, we can obtain an 
anomalously large contribution if there are two distinct reflected rays with nearly the 
same optical path length. If this is the case, it provides the self-consistent justification 
of the approximation. The dominant contribution to the integral will come from 
modes with wavelength of the order of A£; if this is small compared to all other 
length scales in the problem, then the use of geometric optics should be a good 
approximation. 

Normally one would expect Casimir effects to arise from modes whose wavelengths 
are of the order of the length scales defined by the boundaries, typically the distance 
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to the nearest boundary. In this case, one would not expect geometric optics to be a 
good approximation. Nonetheless, Schaden and Spruch || have argued that one can 
often obtain reasonable results from a semiclassical approximation involing a sum over 
periodic classical orbits. Our use and justification of a geometric optics approximation 
is perhaps more akin to that of Hawking || in his derivation of black hole evaporation. 
There the modes which give the dominant contribution to the Hawking radiation have 
very high frequencies when they propagate through a collapsing star, and hence are 
accurately described by geometric optics. 

In this paper, we will examine the case of parabolic mirrors and show that for 
points near the focus, there can be two reflected rays with nearly the same path 
length. Their path lengths differ finitely from that of the incident ray. In this case, 
the dominant contribution to (if 2 ) comes from a single term of the form of T 12 , and 
we can write 

, n, 1 f ,o, COSCo>A£ , . 

^> M 5?/**— — ' (13) 

Note that the interference terms between the incident and the reflected rays give a 
much smaller contribution because the A£ is much larger for these terms. We can 
also now write down expressions for several other quantities of interest. These include 
((p 2 ), where the dot denotes a time derivative, as well as the scalar field energy density 

Pscalar = l^ 2 +\^^\ 2 )-(^). (14) 

In the last step we used the fact that 

l/£°l = |v/£°| (15) 

for plane wave modes. We can also obtain renormalized expectation values for elec- 
tromagnetic field quantities, such as (E 2 ) and (B 2 ), or the electromagnetic energy 
density pem = |((E 2 ) + (B 2 )). Here E and B are the quantized electric and magnetic 
field operators, respectively. The mode functions for these fields are of the form of the 
right-hand-side of Eq. (f|), except with an extra factor of u and a unit polarization 
vector. Thus, when we account for the two polarizations of the electromagnetic field, 
we have 

(E 2 ) = (B 2 ) = Pem = 2(0 2 ) = 2p scalar = — J d 3 ku cos^ A£. (16) 

3 Optics of Parabolic Mirrors 

3.1 Conditions for Multiply Reflected Rays 

A parabolic mirror is illustrated in Fig. 1. The parabola described by 



Figure 1: A parabolic mirror has its focus at the origin. A ray incident at an angle 
9 with respect to the x-axis reflects off the mirror and arrives at a point a distance 
a from the origin at an angle 9'. The line tangent to the point of reflection is at an 
angle a. 
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has its focus at the origin, x = y = 0. Consider a ray incident at angle 9 and reflected 
at angle 9' relative to the x-axis. Further suppose that this ray reaches the x-axis at 
x = a, where a<i). We wish to find the relationship between the angles 9 and 9'. 
First note that 

e = 9'-n + 2a, (18) 

where a is the angle of the tangent to the parabola at the point of intersection. If we 
differentiate Eq. ([17D , we find 



dy_ 

dx 



b 

Vi 



— tan a . 



(19) 



where here yi is the ^/-coordinate of the point of reflection. Note that the reflected 
ray is described by 

y = ttmff (x - a) . (20) 
Combine this relation with Eq. ([[7]) to find 



Vi 



tan 9' 



l±Wsec 2 ^-2(-)tan 2 9> 



(21) 



We expand this expression to first order in a/6 and note that for > 0, we need the 
minus sign before the square root. We then find 



Vi 



tan0' 



Now combine this result with Eqs. ([18]) and (^) to find, to first order in a/6, 

a sin 3 9' sec 9' 
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b(sec9' - 1) 



(22) 



(23) 



First, we note that 9 — > as a — > for fixed 9'. This is the expected result that all 
rays emanating from the focus are reflected into parallel rays. Equation ( ^3l) is plotted 
in Fig. 2. We see that for a ^ 0, there can be two reflected rays for a given incident 
ray. However, one of the reflected rays always corresponds to 9' > n/3. Hence the 
mirror must subtend an angle greater than 7r/3 as measured from the x-axis for this 
to happen. 

Our next task is to compute the difference in path lengths for these two reflected 
rays. Consider first the distance I which a particular ray travels after it first crosses 
the line x = a. This distance can be broken into two segments s\ and S2, as illustrated 
in Fig. 3. If Xj is the x-coordinate of the reflection point, then 



si 



and 



Xi - a) 2 + yt 



Xj — a 



S2 



cos 9 



(24) 
(25) 
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Figure 2: The incident angle 9 as a function of the angle 9' of the reflected ray. Both 
angles are measured in radians. Here it is assumed that the rays arrive near the focus 
(a/b -C 1). For a single incident angle 9, there can be two reflection angles 9'. The 
maximum of this curve occurs at 9' = tt/3. 



We may now use Eqs. ([17]) and (|2T| ) to write 

b 



si 



Similarly, 



Thus, 



sin 2 9' 

b cos 9' 



S2 ~ Xi - a 



sm z 9' 



l-cos0' - (~)sin 2 0' 



S\ + s 2 = b — a(l + cos#') 



(26) 



(27) 



(2* 



If there are two distinct reflected rays with 9' = 9[ and 9' = 9' 2 , respectively, then the 
path length difference is 

A£ = a(cos6[ -cos9' 2 ). (29) 



3.2 Parabola of Revolution 

In Section || we dealt with the rays reflected from a parabola in a plane. There 
are two ways to add on the third spatial dimension. One is to consider a parabolic 
cylinder and the other is to consider a parabola of revolution, the surface formed by 
rotating a parabola about its symmetry axis. In the latter case, one has an azimuthal 



S 



Figure 3: The ray reflects from a point on the mirror with coordinates (xi,yi) and 
then arrives at the point (a, 0). Here S2 denotes the distance traveled by the ray from 
the x = a line to the point of reflection, and si is the distance from that point back 
to x — a. 
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angle < < 2n. Thus Eq. ( |T3"D becomes 

I f r°° I f I 

^-4^I d6 'l ^^o S .A <= -— qb (30) 



Here we have evaluated 



poo poo j 

/ dui uj cos oj Al = lim / du u cosuAl e~ auJ = , (31) 

and then used Eq. fl29l) . 

We can do an analogous calculation for the various quantites given in Eq. (|Tj|) to 
find, for example, 



(E )pr= 2^J (CO8 01- COB 0^ ■ 

Here we have used 

poo poo Q 

/ do; a; 3 cos a; At = lim / <L; uj 3 cos A£ e~ Qa; = — — - . (33) 
Jo a^oJo (A£) 4 K ' 

Note that the integration on 9' in Eqs. (BUI) and (K2|) needs only to include positive 



values of 9'. If one reflects the ray illustrated in Fig. 1 through the x-axis in the 
case of the parabola of revolution, one is going to a ray with the same 9' but with 

(f) —>■ (f) + 71". 

3.3 Parabolic Cylinder 

Another possible geometry in three space dimensions is that of the parabolic cylinder. 
Let the cylinder be parallel to the z-direction. The wavevector k of the light rays 
now has a z-component, k z , so that 



w = ^+A;f, (34) 

where k is the magnitude of the component of k in the xy plane (perpendicular to 
the ^-direction) . If s is a distance traveled in the xy plane, then the actual distance 
traveled is 

a = s-. (35) 

K 

Thus the difference in path lengths for a pair of reflected rays is w Al j k. We can 
modify Eq. ( |T3| ) to give an expression for the mean value of ip 2 near the focus of a 
parabolic cylinder as 

If^^MM. (36) 

87T J J UJ 

Similarly, the mean squared electric field is given by the analog of Eq. (Pf): 

(E 2 ) pc = — f d 3 k u cos(cu 2 A£/k) . (37) 
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The integrations in these two expressions are best done in cylindrical coordinates, 
where d 3 k = ndndO' dk 7 . We can then write 



pc 



8tt 3 
1 

8^ 



de' 



dn k Re e 



ikA£ 



dk 2 



k 2 + k 2 



dO'Re I dKKe^ iKM K (-^iKA£ 



(38) 



where Kq is a modified Bessel function, and in the last step we used Formula 3.364.3 
in Ref. J7|. Next we use Formula 6.624.1 in the same reference to write 



rod 

lim / dx x e~ ax K ((3x) = lim 

a— >fl Jo o— >(5 



a 



P \a 2 -(3 2 V" 2 - P'- 



:ln 



a 



3p 2 



(39) 



We can combine this last result with Eqs. (|30D and ( p8|) to find an expression for 
(^ 2 )pc- However, we need to account for the fact that here, unlike the parabola of 
revolution, the integration on 9' runs over negative values. This can be done by 
introducing a factor of two (corresponding to the contributions of the upper and 
lower halfs of the cylinder) and writing 



Similarly, we can write 
1 



(E 2 ) 



pc 



47T 3 
1 

Arc 3 



If we differentiate Eq. 
that 



(V 2 ) P c = g~ (f 2 )pr ■ 



(40) 



de' 



de' ■ 



dA£ 



dn K 



Im 



d 



dA£ 



dk 7 



: sin 



x - Ik 2 + k 2 



■ 2 i 7,2 \ 
K , 



dKK 2 ei iKM K 



(41) 



with respect to a before taking the limit, we may show 

4 



poo 

lim / dxx 2 e- ax K (px) 

a^{3 Jo 



15/3 

This last identity and Eq. (|i~6"D may be used to show that 

16 



3 " 



pc 



15tt 



(E ) w . 



(42) 



(43) 



Thus the results for the parabolic cylinder are related to those for the parabola of 
revolution by a numerical factor somewhat less than unity. 
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Figure 4: Two flat mirror segments are aligned at angles a\ and a 2 , respectively, and 
subtend angles 0[ and 9' 2 — 0[ from the point of interest. An incident ray has angle 9 
with respect to the x-axis, whereas the reflected ray has angle 9'. 



4 Further Technical Issues 

4.1 Reflected Rays from Line Segments 

In this subsection, we examine the problem of the reflection of rays from a pair 
of attached line segments, as illustrated in Fig. 4. The purpose of this exercise is 
twofold: First, it will lead to the justification of the renormalization prescription 
used in writing down Eq. ([]). Second, it will reveal that integrals such as those in 
Eqs. (|7|) and ( |T6"D should involve an integration over 9', the angle of the reflected 
wave, rather than 9, the angle of the incident wave. 

First consider the case < 9' < 9[, so the ray reflects from the lower segment 
oriented at angle a% relative to the horizontal. Here 9 = 9' + 2a i — n and hence 



Now consider the case where the ray reflects from the upper segment, so 0[ < 0' < 9' 2 
and 9 = 9' + 2a 2 — 7r. Here 



2a x - 7r < 9 < 9[ + 2a\ - 7r . 



(44) 



6[ + 2a 2 - it < 9 < 9' 2 + 2a 2 - n . 



(45) 



Note that the range of 9' 



is AO' = 0' 2 whereas the range of is 



AO = 0' 2 + 2{a 2 - ai) < 0' 2 . 



(46) 
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However, for 9 in the range 9[ + 2a 2 — n < 8 < 8[ + 2a\ — it, there are two reflected 
rays for each incident ray. This is a range of 58 = 2(a 2 — «i), and we have 



58 + A8 = AO' . 



(47) 



Although 8' runs over a larger range than does 8, we can think of this larger 
range as counting the multiple reflected rays that can result from an incident ray 
with a given value of 8. This conclusion will continue to hold if we have more than 
two straight line segments. We can approximate any curve by a sequence of line 
segments. In general, the angle AO' subtended by the curve differs from the range 
of angle of incident rays, A8, and if the curve is convex toward the point of interest, 
A8 < AO'. At first sight, one might think that an emmeration of the independent 
modes should involve an integration over 0. This, however, fails to account for the 
multiple reflected rays, which are correctly counted if we instead integrate on 0' . 

As we vary through its range of —ti < 8 < tt (Note that here 8 increases in the 
clockwise direction.), we have six possibilities: 



In the latter case the incident ray fails to reach the point of interest because it is 
blocked by the mirror. Note, however, that the reflected rays exactly compensate for 
the missing incident rays in the sense that if we add up a weighted sum of the angle 
ranges with reflected rays, it is equal to the range with no rays. This observation 
is the justification for Eq. (Qj. The number of incident plus reflected rays in the 
presence of the boundary is the same as the number of incident rays in its absence. 

One might ask whether it is important also to include the interference terms 
between the multiply refelcted rays. After all, the dominant contribution near the 
focus of a parabolic mirror comes from such an interference term. It is indeed true 
that if one wishes to compute a quantity such as (E 2 ) in the geometry of Fig. 4, 
we would need to include the interference terms. However, one should not expect 
to obtain an anomalously large result, but rather one of order r -4 , where r is the 
distance to the nearest boundary. This follows from the fact that the formula for A£, 
the analog of Eq. ( |2Tf ) will be of the form of a product of r times a dimensionless 
angular dependent function. 

4.2 Evaluation of Singular Integrals 

We have derived expressions, such as Eqs. (|3l|) and (|3"2"|), for renormalized quantities 
near the focus of a parabolic mirror. Recall that we are dealing with a situation 



-7T < 8 < 2«i - 7T 

2«i - 7r < 8 < 9[ + 2a 2 - 7r 
9[ + 2a 2 - vr < 9 < 9[ + 2«i - tt 
9[ + 2ai - tt < 9 < 9' 2 + 2a 2 - tt 
9' 2 + 2a 2 - tt < 9 < 7r - 9' 2 
n-9' 2 <9<n 



incident ray only 



no rays. 



incident ray only 



1 reflected ray 

2 reflected rays 
1 reflected ray 
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where there are two reflected rays for a single incident ray. Here 9[ is the angle of 
one of these rays, and the angle of the other, 9' 2 is understood to be a function of 
9[. However, the integrals in question are singular at the point that 9[ = 9' 2 . The 
singularity may be removed by an integration by parts [§, || jit)] . We rewrite the 
integrand using relations such as 



J_ _lcP_ 

x 2 2 dx 2 



lnx 2 (4£ 



and 



1 1 d 4 , 2 
_ = __ — -lnx 2 . 49 

x 4 12 dx 4 v ; 

Next we perform repeated integrations by parts until we have only an integral with 
a logarithmic singularity in the integrand, plus possible surface terms. Thus, for 
example, 

dx m = jf dxlnx ><?m, ( 50) 

x z 2 J dx z 

and 

dx^ = -± [dxlnx 2 ^, (51) 
x 4 12 J dx A ' v ; 

provided that the function f(x) is regular at x = and the surface terms vanish. 
This procedure is a generalization of the notion of a principal value to cases of higher 
order poles. 

In our case, the integration on 9[ ranges over those values of 9' for which there are 
multiple reflected rays. Within the geometric optics approximation, the integrand 
would seem to drop precipitously to zero at the end point of this interval. If one were 
to go beyond this approximation, the sudden drop would be smeared out over an 
interval corrresponding to about one wavelength. Thus we can think of our integrand 
as being an approximation to a function which, along with its derivatives, vanishes 
smoothly at the endpoints. If so, then we can ignore the surface terms. In any 
case, we will here make the assumption that the surface terms can be ignored. The 
integrand falling smoothly to zero can arise from more than one physical cause. One 
is diffraction effects at the edge of the mirror, as noted above. Another way to enforce 
this behavior is to consider a mirror in which the reflectivity falls smoothly from near 
unity to zero as the edge of the mirror is approached. 

5 Results for Mirrors Slightly Larger than the Crit- 
ical Size 

As we found above (See Fig. 2.), there is a critical size which a parabolic mirror must 
have before we find large vacuum effects near the focus. The critical case is that of a 
mirror which subtends an angle of ir/3 in either direction from the axis of symmetry 
(the x-axis in Fig. 1). In order to evaluate the integrals in Eqs. (|30|) and (|32D , we 
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need to solve Eq. fl23|) for 9' in terms of 9, and then express one root 9' 2 as a function 
of the other, 9[. In general, this is difficult to do in closed form. There is, however, 
one case in which an analytic approximation is possible. This is when the size of the 
mirror is only slightly greater than the critical value. Let the angle subtended by the 
mirror be 7r/3 + £o, where £o "C 1. In this case, we can expand the needed quantities 
in terms of power series. Note that now both roots for 9' will be close to tt/3, so let 
9' = 7r/3 + £ and expand Eq. ( |23"D in powers of £ to find (This and other calculations 
in this section were performed using the computer algebra program MACSYMA.) 

n _3y/3 3^3 2 1 3 3^3 4 1 5 LjV3 6 

U ~— — * +4* 16 4 l80" 4+ "'- 15 j 

Let 6^ = 7r/3 + £i and 9' 2 = vr/3 + £ 2 - Assume a power series expansion for £ 2 in terms 
of £i. Next we equate the right-hand-side of Eq. (|52|) with £ = £i to that with £ = £ 2 
and iteratively solve for the coefficients in the expansion of £ 2 . The result is 

6 - -6 + -y£i - + W £i - ^£i + • • • • (53) 

Our next task is to use this expansion to compute the integrands in Eqs. (P0|) and 
(P^). First rewrite these expressions as 

/pr ~ 87r 2 a 2 7 [cos(f + £ 1 )-cos(§ + £ 2 )] 2 

and 



(V = / ^ 1 r^ | + ei) _ CQs(f , .„„ , (54) 



3 f 1 

(E2)pr = Wi ^ [cos(f + 6) - cos(f + 6)] 4 ■ (55) 

Again we must emphasize that the integrands in these expressions are the approxi- 
mate forms away from the end points of the integrations, but should actually vanish 
at the end points. Next we replace these expressions by the forms obtained by the 
integratios by parts described above, where the surface terms are assumed to vanish. 
After the integrations by parts, we can recognize that the dominant contributions to 
the integrals come from the interval [— £o,£o] and write 



pr 



1 d 2 \ £ 2 1 

—^j^d^lntf^ {[ cos (| + 6 )_ 1 cos (^ + 6 )] 2 ] > {:,b) 



and 



48ir 2 a A J-Ho d£i [ [cos (| + £i) - cos(| + £ 2 J] J 

Note that we have introduced a factor of | to compensate for overcounting of pairs 
of reflected rays. 
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Next we use Eq. (|3|) to write 



[cos(f + a) - cos(| + £ 

and 



A + A 1 £ 1 + A 2 H + ■■■ , (58) 



[cos(f + 6) -cos(f + 



5 + + B 2 £ + £ 3 £ + B£t + ■■■■ (59) 



We see, that to leading order in £o, the dominant contribution to ((f 2 ) comes from 
the coefficient A2, which is given by 



A 2 = — . (60) 
324 K ' 



This leads to our final result 

)pr ~ 648vr 2 a 2 



23 

(A* « -F7^rr£o(l - In^o) + O(£ 2 ln£ ) - (61) 



Similarly, the leading contribution to (E 2 ) pr comes from 

4051 
2 4 3 8 5 



4051 , , 

= 7^ , (62) 



and is 



40^)1 Q 38 x 10~ 3 
(EV « 2 2 3 7 57r2a 4 6>(1 " ln£o) + O(£ 2 ln£ ) « — — A £o(l - ■ (63) 



First we note that the leading contributions to both quantites diverge as a — ► 0, 
that is, as one approaches the focus. This provides the justification of the geometric 
optics approximation. The modes which give the dominant contribution are those 
whose wavelengths are of order a, small enough that geometric optics is valid. Next 
we note that ((p 2 ) pr diverges negatively, but (E 2 ) pr and the energy density for the 
scalar and electromagnetic fields diverge positively. 

The above results apply in the case of a parabola of revolution; for the case of a 
parabolic cylinder we have 

23 

<^—ssv« 1 - 1 »«' (64 » 

and 

< E V - j||^o(l - Ko) « 3 - 18 a X 4 10 ' 3 ;,(l - Ko) • (65) 

Note that all of the results in this section depend upon what is happening in a thin 
band centered on 9' = tt/3. The remainder of the mirror, that for which 9' < 7r/3 — £), 
does not even have to be present. 
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6 Observable Consequences? 



Now we face the question of whether the amplified vacuum fluctuations are actually 
observable. The calculations given above indicate that the energy density and squared 
fields are singular at the focus of a perfectly reflecting parabolic mirror. However, the 
approximation of perfect reflectivity must break down at frequencies higher than the 
plasma frequency of the material in question . So long as the plasma wavelength Xp 
is short compared to the size of the mirror, there is an intermediate regime in which 
geometric optics is valid. We simply must restrict the use of the geometric optics 
results to values of a larger than Xp. 

The quantity which is most easily observable is (E 2 ), as it is linked to the Casimir 
force on an atom or a macroscopic particle. If the atom or particle has a static 
polarizability a, then the interaction energy with a boundary is 

V = ~a(E 2 ). (66) 

Here we are assuming that the modes which give the dominant contribution to (E 2 ) 
have frequencies below that at which a dynamic polarizability must be used. For a 
perfectly conducting parallel plate, 

)platc " 16^ ? ' (6?) 

where z is the distance to the plate. If we insert this expression into Eq. (|66|), then 



the result is the Casimir-Polder potential |LT] for the interaction of an atom in its 
ground state with the plate. It is a good approximation when z is large compared 
to the wavelength associated with the transition between the ground state and the 
first excited state. The 1/z 4 distance dependence of the Casimir-Polder potential was 
experimentally confirmed by Sukenik et al |12| . If we compare Eq. (|67j ) with Eq. (|63|) 



or Eq. (|6q), we see that the mean squared electric field near the focus of a parabolic 
mirror is only slightly less than that at the same distance from a flat plate. Given that 
the latter has actually been observed, it is possible that the inhanced fluctuations 
near the focus are also observable by techniques similar to those by Sukenik et al. 

The basic method used in the Sukenik et al experiment is to look for the effects 
of the deflection of a beam of atoms as it passes near a pair of plates. We can give 
a general estimate of the size of this type of deflection which applies whenever there 
is a mean squared electric field which varies as the inverse fourth power of a length 
scale. Let 

(E 2 ) = 4 , (68) 
or 

where a is the length scale and A is a dimensionless constant. We assume that an 
atom has an interaction of the form of Eq. (0). The resulting force, F = —dV/da, 
will cause a deflection Aa in the atom's position in a time t, where 

(^)(i)(^)(^) 6 (i^) 2 - w 
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Here m^ a = 3.8 x 10~ 23 gm and a^ a — 3.0 x 10~ 22 cm 3 denote the mass and polar- 
izability of the sodium atom, respectively. (Note that polarizability in the Lorentz- 
Heaviside which we use is An times that in Gaussian units.) If t is of order 10~ 3 s 
(the time needed for an atom with a kinetic energy of order 300K to travel a few 
centimeters), and z is of order 1/im, the fractional deflection is significant. Recall 
that in our case 

| 9.38 x lCT 3 £o(l — ln£o) , parabola of revolution , . 

| 3.18 x 10~ 3 £ (1 — hi£o) , parabolic cylinder. ^ ' 

Thus it may be possible to observe the force on atoms near the focus. 

Another possible way to observe this force might be to levitate the atoms in the 
Earth's gravitational field. (A rather different form of levitation by Casimir forces 
was proposed in Ref. @.) If one equates the force on atom at a distance a from the 
focus to its weight, the result can be expressed as 

Given that this formula applies for a > Xp and that Ap ~ 0.1/im for many metals, it 
seems possible that levitation near the focus is possible. Of course, atoms will only 
be trapped if their temperature is sufficiently low. The required temperature can be 
estimated by setting the thermal energy |fcT equal to the magnitude of the potential 
energy V. The result is 

Thus for a of the order of a few times 0.1/i.m, the required temperature is larger than 
the temperatures of the order of 10~ 7 K which have already been achieved for laser 
cooled atoms |L3|, |i4| . 

Another possibility might be the use of atom interferometry. Atoms traveling for 
a time t parallel to and near the focus of a parabolic cylinder will acquire a phase 
shift of 

^=i a <EV = 0.14 (£, (l^) 4 ( T5 L ; ) 6 (1 -!»&). (73) 

If it is possible to localize the atoms to within a few fim of the focus, then the 
accumulated phase shift for reasonable flight times would seem to be within the 
currently attainable sensitivities of the order of 10~ 4 radians fT3|. 



7 Discussion and Conclusions 

In this paper we have argued that a parabolic mirror is capable of focusing the vacuum 
modes of the quantized electromagnetic field and creating large physical effects near 



18 



the mirror's focus. Just as the mirror can focus a beam of light, it can focus something 
even in the absence of incoming light. This might be dubbed "focusing a beam of 
dark" ||15|| . The manifestation of this focusing is a growth in the energy density and 
mean squared electric field as the focus is approached. In the idealized case of a 
perfectly reflecting mirror, these quantities diverge as the inverse fourth power of the 
distance from the focus. For a real mirror, the growth is expected to saturate at 
distances of the order of the plasma wavelength of the mirror. 

The most readily observable consequence of the focused vacuum fluctuations is 
enhanced Casimir forces on atoms or other particles near the focus. The sign of the 
force is such as to draw particles into the vicinity of the focus. Estimates given in 
the previous section indicate that the magnitude of this effect may be large enough 
to be observable. 

The calculations presented in this paper were based on the geometric optics ap- 
proximation in which only short wavelenth modes are considered. The justification of 
this approximation is self-consistency: the large effects near the focus can only come 
from the short wavelength modes for which the approximation should be a good one. 
Nonetheless, in future work it will be of interest to go beyond the geometric optics 
approximation. This should allow one to check the validity of the assumption made 
in Sect. |4.2| that the surface terms can be ignored. 

In order to simplify the calculations, we made two restrictions on the geometry. 
The first is that we have assumed that the point at which the mean squared field 
quantites are measured lies on the symmetry axis of the parabola (the x-axis). The 
second is that the mirror be only slightly larger than the critical angle of 7r/3 at which 
vacuum focusing begins. (This is the assumption that £ <C 1, made in Sect. ||.) It is 
of interest to remove both of these restriction, which we hope to do in a future work. 
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